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Abstract 

The aim of this paper is to extend the approximate quasi-interpolation on a uniform grid 
by dilated shifts of a smooth and rapidly decaying function on a uniform grid to scattered 
data quasi-interpolation. It is shown that high order approximation of smooth functions 
up to some prescribed accuracy is possible, if the basis functions, which are centered at the 
scattered nodes, are multiplied by suitable polynomials such that their sum is an approximate 
partition of unity. For Gaussian functions we propose a method to construct the approximate 
partition of unity and describe the application of the new quasi-interpolation approach to 
the cubature of multi-dimensional integral operators. 

1 Introduction 



The approximation of multivariate functions from scattered data is an important theme in 
numerical mathematics. One of the methods to attack this problem is quasi-interpolation. One 
takes values M(xj) of a function n on a set of nodes {xj} and constructs an approximant of u by 
linear combinations 

^M(xj)r/j(x), 

where ?7j(x) is a set of basis functions. Using quasi-interpolation there is no need to solve large 
algebraic systems. The approximation properties of quasi-interpolants in the case that xj are 
the nodes of a uniform grid are well-understood. For example, the quasi-interpolant 

E-(^J)^(^) (1-1) 

can be studied via the theory of principal shift-invariant spaces, which has been developed in 
several articles by de Boor, DeVore and Ron (see e.g. [2], 0)- Here ip is supposed to be 
a compactly supported or rapidly decaying function. Based on the Strang-Fix condition for 
which is equivalent to polynomial reproduction, convergence and approximation orders for 
several classes of basis functions were obtained (see also Schaback/Wu '201, Jetter/Zhou [71). 
Scattered data quasi-interpolation by functions, which reproduce polynomials, has been studied 
by Buhmann, Dyn, Levin in and Dyn, Ron in 0j (see also |25 for further references). 



In order to extend the quasi-interpolation Hl.lj) to general classes of approximating functions, 
another concept of approximation procedures, called Approximate Approximations, was proposed 
in 1^ and |lL)j . These procedures have the common feature, that they are accurate without being 
convergent in a rigorous sense. Consider, for example, the quasi-interpolant on the uniform grid 

M n(x) = D-^/' um , (1.2) 

where t] is sufficiently smooth and of rapid decay, h and D are two positive parameters. It was 
shown that if J^rj — 1 has a zero of order N at the origin (J^rj denotes the Fourier transform of 
rj), then M u approximates u pointwise 

\Mu{x) - u{x)\ ^ CN,r, {hVD)^ sup \Vnu\ + e \Vn-Mx)\ (1.3) 

with a constant cn^t^ not depending on u, h, and D, and e can be made arbitrarily small if D 
is sufficiently large (see JHl ) !14j ) . In general, there is no convergence of the approximate quasi- 
interpolant A4u{x) to ?i(x) as /i — > 0. However, one can fix D such that up to any prescribed 
accuracy Aiu approximates u with order 0{h^). The lack of convergence as /i ^ 0, which is 
even not perceptible in numerical computations for appropriately chosen D, is compensated by 
a greater flexibility in the choice of approximating functions 77. In applications, this flexibility 
enables one to obtain simple and accurate formulae for values of various integral and pseudo- 
differential operators of mathematical physics (see ^21; ^Hj; QZ] and the review paper 21 J and 
to develop explicit semi-analytic time-marching algorithms for initial boundary value problems 
for linear and non linear evolution equations [Hj). 

The approximate quasi-interpolation approach was extended to nonuniform grids up to now 
in two directions. The case that the set of nodes is a smooth image of a uniform grid have 
been studied in . It was shown that formulae similar to H1.2|l preserve the basic properties of 
approximate quasi-interpolation. A similar result for quasi-interpolation on piecewise uniform 
grids was obtained in [Hj. 

It is the purpose of the present paper to generalize the method of approximate quasi- 
interpolation to functions with values given on a rather general grid. We start with a simple 
quasi-interpolant for a set of nodes close to a uniform grid of size h in the sense, that for some 
positive constant k and any j G there exists at least one node xj with |xj — < nh. Then 
under some additional assumption on the nodes we construct a quasi-interpolant with gridded 
centers 

Here Aj are linear functionals of the data at a finite number of nodes around xj. It can be shown 
that estimate remains true for Mn under the same assumptions on the function 77. 

In order to treat more general distributions of the nodes xj we modify the approximating 
functions. More precisely, we consider approximations of the form 

Mn(x) = ^^(xj)Pj(x)r,(^i^) (1.4) 
j ^ 

with some polynomials Vy We show that one can achieve the approximation of u with arbi- 
trary order AT up to a small saturation error, as long as an "approximate partition of unity" 

|7'j(x)}r7^ — ^- — I with other polynomials exists. Here we mean that for any e > one 
can find polynomials such that 



< e. 



Then one can choose the polynomials Vj in H1.4|) such that 

\Mu{^) - n(x)| ^ C sup /if \\Vnu\\l^ + e |n(x)| . 
j 

This estimate is valid as long as 

j 

and r] is sufficiently smooth and of rapid decay, but is not subjected to additional requirements 
as the Strang-Fix condition. Moreover, we propose a method to construct the polynomials such 
that the series 

j 

approximates the constant function 1 up to an arbitrary prescribed accuracy. This method 
does not require solving a large system of linear equations. Instead, in order to obtain the 
local representation of the partition of unity, one has to solve a small number of approximation 
problems, which are reduced to linear systems of moderate size. 

By a suitable choice of r] it is possible to obtain explicit semi-analytic or other efficient 
approximation formulae for multi-dimensional integral and pseudo-differential operators which 
are based on the quasi- interpolant H1.4|) . So the cubature of those integrals, which is one of the 
applications of the approximate quasi-interpolation on uniform grids, can be carried over to the 
case when the integral operators are applied to functions given at scattered nodes. 

We give a simple example of formula ()1.4|) . Let {xi} be a sequence of points on M such that 
< Xj+i — Xi ^ 1. Consider a sequence of functions Cj on R supported by a fixed neighborhood 
of the origin. Suppose that the sequence {Cjix — Xj)} forms an approximate partition of unity 
on M, 

\l-^Cjix-Xj)\ <e . 
j 

One can easily see that the quasi-interpolant 
satisfies 

\Mhu{x) - u{x)\ ^ch^ II^^"IIloo{M) +£\u{x)\ , 

where the constant c depends on the functions 

The outline of the paper is as follows. In Section[21we consider an extension of the approximate 
quasi-interpolation to scattered nodes close to a uniform grid. We construct the quasi-interpolant 
Mu with gridded centers and coefficients depending on scattered data and obtain approximation 
estimates. Further the results of some numerical experiments are presented which confirm 
the predicted approximation orders. In Section we show that an approximate partition of 
unity can be obtained from a given system of rapidly decaying approximating functions if these 
functions are multiplied by polynomials. Using the approximate partition of unity, one can 
construct approximate quasi-interpolants of high order approximation rate up to some prescribed 
saturation error. This is the topic of Section [l] Section [5] contains an application to the cubature 
of convolution integral operators. A construction of the approximate partition of unity for the 
case of Gaussians and some numerical examples are given in Section |HJ 



2 Quasi-interpolants with gridded centers 

Here we give a simple extension of the quasi-interpolation operator on uniform grids H1.2() to a 
quasi- interpolant, which uses the values u(xj) on a set of scattered nodes X = {xj} C M" if it is 
close to a uniform grid. Precisely we suppose 

Condition 2.1 There exist h > and ki > such that for any j S the ball B(hj,hKi) 
centered at h} with radius hKi contains nodes of X. 



{N -1 + n) 
n\{N-l) 

called the star of xj and denoted by st (xj) if the Vandermonde matrix 



Definition 2.1 Let x; G X. The collection of uin = — — — 1 nodes Xk G X will be 

■> nU A/ — 1 ll 



= l"l = l'-'^-l' (2-1) 

is not singular. The union of the node xj and its star is denoted by ST (xj) = xj U st (xj). 

Condition 2.2 Assume Condition 12.11 and denote by xj G X the node closest to hj. There 
exists K2 > such that for any j G 

(a) st (xj) C -B(xj, h K2) with | det Vj,h\ ^ c > 0; 

(b) U ST(5j) = X. 

2.1 Error estimate 

To formulate our first result we denote by {^^^^i' -^k S st (xj), \cx\ = 1, ...,N — 1, the elements 
of the inverse matrix of V^^h and define the functional 

N~l ~ N-l ~ 

|ql|=1 Xi56st{xj) XkSst(xj) |ct|=l 

Theorem 2.1 Suppose that for some K > n and the smallest integer uq > n/2 the function 
r/(x), X G M", satisfies the conditions (1 + |x|)^|(9^7/(x)| < for all < |/3| < uq, and 
d"^{J^rj — 1)(0) = 0, < |q;| < A'". // the set of nodes X satisfies Conditions 12.11 and 12. 2L then 
for any e > there exists D such that the quasi-interpolant 



lu 



(x) = Z?-«/2^A.(,)^(^) (2.2) 



approximates any u G (M") with 

N-l 

|Mu(x) - n(x)| ^ CNr^D /i^supIVatuI + e V |Vfctt(x)|(\/7^/i)^ , (2.3) 
where CN,r),D does not depend on u and h. 

Proof. For given u G VF^(M") we consider the quasi-interpolant (|1.2)) on the uniform grid {/ij} 

'y.-hy 



with h given by Condition l2.1l It was proved in jl6j that under the decay and moment conditions 
on r], formulated in the statement of the theorem, A4 u can be represented as 

A^n(x) = n(x)+ g (^)'"'^ 9'^^,(./^k) e^C^-) + f/.(x) 

|a|=0 ' kGZ"\{0} 

with a function Un bounded by 

\Un{^)\ <Ci {^Dh)^ sViv\V Nu\ 

and a constant Ci depending only on rj. Moreover, the sequences {d'^Tr}{\/lD •)} G and 

^ |5"J^?7(/Dk)| ^0 as D ^ oo . 
keZ"\{o} 

Hence, we can find D such that ^Au satisfies the inequality 

N-l 

|7Wu(x) -tx(x)| ^ \UN\+eY,\^ku{-K)\{^hf . 

k=0 

It remains to estimate \Mu — Mn|. Recall the Taylor expansion of u around y G 



W= E ^^(x-yr + i?7v(y,x) (2.4) 

|a|=0 



with the remainder satisfying 

|i?Ar(y,x)| ^ CAr|x - y|^ sup \Vnu\. (2.5) 

B(y,|x-y|) 

For j G we choose xj G X and use (|2.4() with y = xj . We split 
7Wu(x) = M(i)n(x) + i?jv(Xj, /ij) 

with 

Af(')„w = «-»/'^X: (2.6) 

jez" |a|=o ■ ^ 
Because of | /ij — xj | < ki /i for any j we derive from ()2.5() 

|M(i)n(x) - ^^^(x)! ^ CAT (ki/i)^D-"/2 Y Lf^LZ^\ gup |v^u| . (2.7) 



B{x,hKi) 



The next step is to approximate 9"u(xj), 1 ^ |q:| < A^, by a linear combination of u(xk 

;|«| 



Xk G st (xj). Let {aii'*}i<g|Q,|<jv-i be the unique solution of the linear system with ttin unknowns 



(j) 

E ^j^" " ~ ^^^i^ ' e St (Xj) . (2. 

I«l=i 



From and (|^ follows that 



|«|=i 



a! 



i?Ar(xj,Xk) . 



By Condition 12. 2f a) the norms of are bounded uniformly in j, this leads together with (|2.5I) 
to the inequality 



■J^' ^ C2/i^-l"l sup |V^u|. 

_B(Xj,/lK2) 

(j) 



Hence, if we replace the derivatives (?"u(xj) in (|2.6() by aa , then we get the sum 



l«|=i 



which in view of 



«-=7^ E ^SkKxk)-n(5j)) 



XkGst (xj) 



coincides with H2.2() . Moreover, 



(x) - M(i)tx(x)| ^ Cs/i^ ^ 

|«|=i 



X - /ij 



h^D 



sup IVatuI . 

B(x,hK2) 



x-J 



Now the inequality 

supL>-"/2 yu 

for all D > Do > implies that (fT7)) and (fTTU)) lead to 

|>/u(x) - Mn(x)| < Q/i^ sup |Vjvu| , 

B(x,/lK2) 



(2.9) 



(2.10) 



which proves (|2.3jl . 



2.2 Numerical Experiments with Quasi-interpolants 

The behavior of the quasi-interpolant Mu was tested by one- and two-dimensional experi- 
ments. In all cases the scattered grid is chosen such that any ball B{h],h/2), j € Z", con- 
tains one randomly chosen node xj, thus xj = xj. All the computations were carried out with 
MATHEMATICA®. 

The following figures show the graph of M.u — u for different smooth functions u using basis 
functions for second (Fig. ^ and fourth (Fig. |2I) order of approximation with h = 1/32 (dashed 
line) and /i = 1/64 (solid line). 

In Fig. 13] the difference Mn — n is plotted for the function u(x) = (1 + |x|2)-i, X G R"^. 

In formula ()2.2() we use 7/(x) = 7r~^e~l''l , D = 2, for which N = 2, and therefore the star 
st(xj) contains 2 nodes. We have chosen st{xj^jr^) = {xj^+ijj, Xj^ jj+i}. In Tabled we give 
some computed values of Mm(0) — u(0) for D = 2 and D = 4 with different h, which confirm 
/i^-convergence of the two-dimensional quasi-interpolant. 



Figure 1: The graphs of Mu{x) — u{x) with D = 2, N = 2, st (xj) = {xj+i}, when u{x) = (on the 
left) and u{x) = 1/(1 + x^). Dashed and sohd lines correspond io h = 1/32 and h = 1/64. 
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Figure 2: The graphs of MM(a;) — u{x) with £) = 4, = 4, st {xj) = {xj-2,Xj-i,Xj+i}, when u{x) = x^ 
(on the left) and u{x) — 1/(1 + x^). Dashed and solid lines correspond to h = 1/32 and h = 1/64. 



3 Approximate partition of unity 

In the following two sections we consider irregularly distributed nodes. First we show that an 
approximate partition of unity can be obtained from a given system of approximating functions 
centered at the scattered nodes if these functions are multiplied by polynomials. We are mainly 
interested in rapidly decaying basis functions which are supported on the whole space. But we 
start with the simpler case of compactly supported basis functions. 



3.1 Basis functions with compact support 

Lemma 3.1 Let {i?(xj, /ij)}j^o be an open locally finite covering ofW^ by balls centered in xj 
and radii hj . Suppose that the multiplicity of this covering does not exceed a positive constant 
Hn and that there are positive constants ci and C2 satisfying 



r 



(3.1) 



provided the balls B(xj,hj) and B{x^, hai) have common points. Furthermore, let {r/j} be a 
bounded sequence of continuous functions on R" such that supp?7j C i?(xj, /ij). We assume that 
the functions 9 y — ^ ??j (^j y ) are continuous uniformly with respect to j and 



«W '-^Vii^) ^ c on 



(3.2) 



Figure 3: The graph oiMu-u with D = 2, N = 2, h = 1/128, when u(x) = 1/(1 + [xp), x e R^. 
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Table 1: Values of Mu(0) - u{0) 



where c is a positive constant. Then for any e > there exists a sequence of polynomials {V^} 
with the following properties: 

(i) the degrees of all Vj are bounded (they depend on the least majorant of the continuity 
modulae of rjj and the constants e, c, c\, C2, /x^j; 



(ii) there is such a constant cq that |7^j| < Cq on i?(xj,/ij); 

(iii) the function 



(3.3) 



satisfies 



|e(x)-l|<£ /ora//x€M", 



(3.4) 



Proof. Since the functions B{xj, 1) 3 y ^ s{hj y) are continuous uniformly with respect to j, 
for an arbitrary positive 6 there exist polynomials T-'j subject to 



.(x) 



< S on B(xj, /ij) , 



and the degree of Pj, degVj, is independent of j. Letting 5 = e {fJ-nWvWLoo) ^ '^^ obtain 



s(x) 



< 



Then 



<e. 



since at most /i„ terms of this sum are different from zero. But 

j j j j 

which proves (|3.4() . 



(3.5) 
(3.6) 



Remark 3.1 Let the functions {r/j} in Lemma f3.1l satisfv the additional hypothesis ryj G C'^(M") 
Then one can find a sequence of polynomials {'Pj} of degrees Lj such that 



sup 



^j(x) 



.(x) 



-j iJ(Xj,ftj) 



^C(A:)-^ sup |Vfcs(x) 



(see, e.g., This shows that it suffices to take polynomials "Pj with deg'Pj > c{k)£ in 

order to achieve the error e in ()3.6() . 

3.2 Basis functions with noncompact support 

Here we consider approximating functions supported on the whole M". We suppose that the 
functions r/j are scaled translates 

of a sufficiently smooth function r/ with rapid decay. 

Lemma 3.2 For any e > there exists and polynomials {"Pj} of degree degPj < such 
that the function Q defined by (|3.3() satisfies 1)3. 4|) under the following assumptions on rj, the 
nodes {xj} and the scaling parameters {hj}: 



1. There exists K > such that 

:=||^(l + /,-i|.-Xj|)- 



K 



< oo . 



(3.7) 



2. There exists p > such that 

,K „2|.|2 



(l + |.|)^e^''l-l\ , (l + |.|)^Vr? 



< Cp < oo . 



3. There exists C > 1 such that all indices j, m 



< C . 



(3.8) 



(3.9) 



4. is valid. 



Proof. Prom p.7|) and (|3.8j) the sum s(x) = jyj (x) converges absolutely for any x to a positive, 
smooth and bounded function s. Suppose that we have shown that for any e > and all indices 
j there exist polynomials Vj such that 



1 



.(x) 



CK ^ ■' ■' ^ 



-K 



{ck is defined in ((213)) and degPj < L^. Then 

sup^ Lj(x) m(x 



.(x) 



and as in the proof of Lemma 13.11 we conclude 

sup K]]r?j(x)-Pj(x) - 1 



< e. 



(3.10) 



To establish ()3.10() we use a result on weighted polynomial approximation from [3], which will 

2 I 1 2 

be stated in a simplified form. Introduce the weight Wp{:x.) = , p > 0, and consider the 

best weighted polynomial approximation of a function g given on 



EN{g)wp 



,oo := inf sup \wp{x){g{x) - P(x))| 



where Hn denotes the set of polynomials, which are of degree at most N in each variable 
xi,...,Xn. Then for g G W^{W) 



(3.11) 



Let us fix an index j and make the change of variables y = /ij ^(x — xj). Then H3.10|) is proved 
if we show that there exists a polynomial "Pj such that for all y G 



«(y) 



<-(i + |y|)-" 

CK 



(3.12) 



with s(y) = s(/ijy + xj). Since s ^ £ W^{W^) according to (|3.11|) we can find a polynomial "Pj 
satisfying 



sup 



< 



' s(y) CpCK 

with the constant Cp in the decay condition 1)3. 8() . Now 1)3. 12() follows immediately from 

\v{y)\{i + \y\f <cpe-p'\yr 

From (|3.11|1 we see that deg Vj depends on the norm of the gradient 



sup 



V 



s(y) 



h] sup 



V 



.(x) 



1 \ - hi 

< sup , , > - — 



X — Xr 



which is bounded uniformly in j in view of 1)3.2(1 , (|3.8|) , and (|3.9() . 



4 Quasi-interpolants of a general form 



In this section we study the approximation of functions u G W^(M") by the quasi-interpolant 
H1.4|) . We win show that within the class of generating functions of the form polynomial times 
compactly supported or rapidly decaying generating function it suffices to have an approximate 
partition of unity in order to construct approximate quasi-interpolants of high order accuracy 
up to some prescribed saturation error. 

Let us assume the following hypothesis concerning the grid {xj}: 

Condition 4.1 For any xj there exists a ball B(xj,hj) which contains m^r nodes Xk € st (xj) 
with 

|de.V5,,| = |de,{(2i-i!i)°};:jj>.. (4.1) 



(see Definition l2.1|) . with c > not depending on xj. 



4.1 Compactly supported basis functions 

Theorem 4.1 Suppose that the function system {rj^} satisfies the conditions of Lemnia \'6.1l let 
u G W^{W^) and e > arbitrary. There exist polynomials Vj,k, independent on u, whose degrees 
are uniformly hounded, such that the quasi-interpolant 

Mu{x) = Y,n{x^) Yl ^j-k(x)^j(x) (4.2) 

k ST(xj)9Xk 

satisfies the estimate 

\Mu{x) -u{x)\ s^Ch^ sup \Vnu\ +£\u{x)\, (4.3) 

where x^ is an arbitrary node and x is any point of the ball i?(xin, /im). By A we denote a 
constant greater than 1 which depends on c\ and C2 in 1)3. The constant C does not depend 
on km, m and e. 

Proof. For given e we choose polynomials 'Pj(x) such that the function satisfies 

|e(x) - 1| < e for all x G R" , 
and introduce the auxiliary quasi-interpolant 

M(i)^(x) = E ( E -^^^ - Xj)")^j(x)r?j(x) . (4.4) 
j |a|=0 

Using the Taylor expansion (|2.4j) with y = xj we write M(^)ti(x) as 

M«^x(x) = n(x)e(x) - ^i?^(xj,x)Pj(x)77j(x) , 

j 

which gives 

|M«t.(x) - n(x)| ^ l^iv(xj,x)Pj(x)7?j(x)| + |n(x)| |e(x) - 1| . 
j 

This, together with the estimate for the remainder ()2.5() . shows that for x G i?(xjn, km) 

|M(i)n(x) -^(x)! ^ C7i/i^ sup | Vjv u| + e|n(x)| , (4.5) 

B(Xin,A hm) 



where the ball i?(xm,A/ini) contains all balls -B(xj, hj) such that -B(xj, hj) and -B(xm, /im) 
intersect. 

Similar to the proof of Theorem 12.11 we approximate in M^^^u the values of the derivatives 
9"ti(xj) by a linear combination of M(xk), where x^ S st (xj). The solution of 

XI -^i^^- ^j)" = ^(^k) - ^^(Xj) , Xk e st (xj) , 
|«|=1 



is given by 



"j Xkest(xj) 



where {^^^^i elements of the inverse of V^^hy Replacing the derivatives {9"u(xj)} in 

by {a£^} gives the quasi-interpolant 

N-l 

M.(x) = ;X{.(x,)(i- Ee(^)" 

j Xkest(xj) |«|=i J 

Xkest(xj) l"l=i 
= E E ^(xk)7^j,k(x)r/j(x) 

j XkeST(xj) 

which can be rewritten as the quasi-interpolant (|4.2)1 . By (|2.4j) we obtain again 

N-i M _ 
E ^(a^^-5-u(xj))(^^)" = i?^(xj,Xk), 
|«|=i ■ J 

hence the boundedness of ||V.~,^|| from Condition 14.11 and the estimate of the remainder (|2.5|1 
imply 

|aH^ -a"n(xj)| < QlCa/if"'"' sup \Vnu\. 

iJ(Xj,ftj) 



Therefore we obtain the inequality 



N-l I I 



|Mu(x) -M«'u(x)K C79 V/if sup Nnu\S2\—-^ |Pj(x)?7j(x) 



and, for any x G B{x^, h^), 



|Mn(x) -M(i)u(x)K Cg/i^ sup |V7vm|. 

B(Xi„,A/lm) 



This inequality and 1)4. 5|) lead to 1)4. 



4.2 Quasi-interpolants with noncompactly supported basis functions 
Theorem 4.2 Suppose that additionally to the conditions of Lemma 13.21 the inequality 



E(i + Vi--xji) 



N-K 



< OO 



(4.6) 



is fulfilled, let u G and e > arbitrary. There exist polynomials V^^]i_, independent on 

u, whose degrees are uniformly hounded, such that the quasi-interpolant 



k ST(xj)9Xk 



X — X; \ /X — X; 



h^ 



hi 



satisfies the estimate 



I Mti (x) — ii(x)| ^ C sup /ij^ llVArnlli^ +e|ii(x 



(4.7) 



The constant C does not depend on u and e. 

Proof. Analogously to (|4.4|) we introduce the quasi-interpolant 



M<'>„(x)=j:(5:^^(.-.ir)T'i 



X — X;\ /X — X 



j |a|=0 



and obtain the estimate 



|MWu(x) -^(x)! ^ ^|i?iv(Xj,x)Pj 



X — X; \ /X — X 



From (|3.12|) we have 



X — X; \ /X — X; 



1 

< - 

c 



h- 



X — X; 



+ |^x(x)||e(x)-i| 



+— (1 + 

CK 



|x-Xjh-i^ 



with the lower bound c of s(x) (see (|3.2|) 1. Together with (|3.8|1 and (|2.5j) this provides 



-RAf(xj,x)'Pj 



X — X; \ /X — X; 



"•"^ H h- 

< CNhf\\VNu\\L„ 



X — X; 



— e 

c 



CK 



-K 



resulting in 

|mWu(x) - n(x)| < |u(x)| |G(x) - 1| + cn\\Vn u\\l, 
^||e-f'W'|xr|L y/iffl + 



J 



N-K 



Now we can proceed as in the proof of Theorem 14.11 

Remark 4.1 Let for fixed x the parameter Kx be chosen such that 



E 



Vlx-xjiV/^l 



X — X; 



|Xj-x|>Kx 



N 



1 + 



X — Xi 



-K 



< £. 



Then the estimate (|4.8|) can be sharpened to 



\Mu{x) — u{x.)\ ^ C max h: sup IVatuI + e (|^i(x)| + || Vat 

|Xj-x|<Kx B(x,Kx) 



5 Application to the computation of integral operators 

Here we discuss a direct application of the quasi-interpolation formula 1)4 .71) for the important 
example r/(x) = e^l''' . Suppose that the density of the integral operator with radial kernel 

/Cn(x) = 1 5(|x-y|)n(y)dy (5.1) 

is approximated by the quasi-interpolant 

Mn(x) = j; ^(xk)^J,k(^)e-l— ^l^/^r (5.2) 

j Xk6ST(xj) J 

Using the following lemma it is easy to derive cubature formulae for (|5.1|1 . 

L L 

Lemma 5.1 For any P(x) = c^x^ one can write P(x)e~'''' = S /^{dx) e~^^^ with 

|/3|=0 |/3|=0 
the polynomial 5/3 (t) being defined by 

^Mt)^(^)""«.(^), (-) 

where Hf^ denotes the Hermite polynomial of n variables Hf^{t) = el*l^(— c?t)^e~l*l^. 
Proof. We are looking for the polynomial 5^(t) defined by the relation 

x/3e-|xp = 5/3(5x)e-l^l', X e M". (5.4) 

Since 

^(5/3(ax)e-l'^l')(A) = W2e-"'l^l'5/3(27riA) 

and 



^(x^e-W^)(A) = W2(- A)V-^I^I 



2| \|2 



we obtain (|5.:-{j) . 



In view of Lemma 15.11 we can write Pj,k(x)e-I^l' = 7j,k(5x) e"!''!' with some polynomials 
7j^k(x). Then 1)5.21) can be rewritten as 

M^(x) = ^(xk) 7j,k(-/ij a..) e-l'^-^jl'/^' . 

j XkeST{xj) 

The cubature formula for the integral Ku is obtained by replacing u by its quasi-interpolant Mu 
£n(x) = /CM^x(x) = Y, E ^(xk)'?],k(-^j5xj)/ij^ y 5(/ij|z|)e-l^+tjl'dz, (5.5) 

j XkeST(xj) j5„ 

where tj = (x — Xj)//ij. By introducing spherical coordinates in we obtain 



oo 



y5(^j|z|)e-l^+*jl'dz = e-l*jl'y" ^>""^5(/ij^) e'^'ci ^ y" 



5^ 



where ^ is the unit sphere in W^. The integral over ^ can be represented by means of 
the modified Bessel functions of the first kind In in the fohowing way 

TT 

/9 (n— 1)/2 r 

Si-i 

(see m p. 154] and |2S1 p.79]). If we denote by 

oo 

£(r) = 2 W2 ri-"/2e--' ^ f,"/^ e-^'5(/ij q) hn-mi^gr) dg , 



then (|5.5j) leads to the following cubature formula for the integral Ku 



j XkeST(xj) J 

6 Construction of the 0-function with Gaussians 

In this section we propose a method to construct the approximate partition of unity for the 
basis functions 

?7j(x) = (vr iJ) ' e ' J' ' j 
if the set of nodes {xj} satisfy Condition 12.11 piecewise with different grid sizes hj. 

6.1 Scattered nodes close to a piecewise uniform grid 

Let us explain the assumption on the nodes: Suppose that a subset of nodes xj G Jq satisfies 
Condition 12.11 with h = hi. The remaining nodes Xk G X \ Jq lie in a bounded domain 17 C M" 
and satisfy Condition 12.11 with /i = /12 = Hhi for some small H. To keep good local properties 
of quasi-interpolants one wants to approximate the data at these nodes by functions of the form 
polynomial times e"'''"''"'' /^^a^^ whereas outside Q quasi-interpolants with functions of the form 
polynomial times e~l^~-^jl Z'*!^ should be used. 

Our aim is, to develop a simple method to construct polynomials Vj such that 

e(x) = (vrz?)-«/2( y: Hlrmy^'^'^^^'''' + ^ ^'^(ir^i)^"'™''^'''') ^'-'^ 

is almost the constant function 1. Here J2 denotes the set of nodes x^ € and Ji = {xj} \ J2 
the remaining nodes. 

First we derive a piecewise uniform grid on which is associated to the splitting of the set 
of scattered nodes into Ji and J2. We start with Poisson's summation formula for Gaussians 

which shows that 



<Cie--'^ 



with some constant Ci depending only on the space dimension. 



Thus for any e > there exists D > such that the function system {e~^^~^^"^^^ 
forms an approximate partition of unity with accuracy e. We can represent any of these functions 
very accurately by a hnear combination of dilated Gaussians due to the equation (see [TS]) 



E 

(6.2) 

-e-l'^lV^i ^ ^2TTi{Di-h^D)i^,k)/hDi ^-7T^D{Di-h^D)\k\^/Di^ 

keZ"\{o} 

which is valid for any Di > /i^Z) > 0. Applied to our setting with h = h2 and Di = h\D we 
obtain the approximate refinement relation 

^-\.\VhlD _ Y ake-l'^-'^^'^l'/'^^^l < C2e-l-l'/^?^e-'^(i-^') (6.3) 
(because by assumption /i2 = Hhi) with the coefficients 

Again, the constant C2 depends only on the space dimension. Define by 5* G the minimal 
index set such that 

Ok < e ^ ^ 

keZ"\s 

Then it is clear from 1)6. 3|) that for any disjoint Zi and Z2 with ZiU Z2 = 

meZi meZ2 keS 

(6.4) 

Condition 6.1 Denote Z2 = {m G : /iim + /12k G for all k G S"}. The constant ki of 
Condition 12. II and the domain are such that for all nodes Xk G 17, i.e. the nodes belonging to 
J2, one can find m G Z2, k G S" with |xk — him — /i2k| < Ki/i2. 

Setting Zi = Z" \ Z2 we connect the index sets Zi , Z2 with the splitting of the scattered nodes 
into Ji, J2. By this way we construct an approximate partition of unity using Gaussians with 
the "large" scaling factor hi centered at the uniform grid Gi := {/iimlme^i outside and using 
Gaussians with scaling factor /12 and the centers G2 ■= {hiui + /i2k}meZ2,keS ^• 

It is obvious, that the above definition of piecewise quasi-uniformly distributed scattered 
nodes and the construction of an associated approximate partition of unity on piecewise uniform 
grids can be extended to finitely many scaling factors hj. Since there will be no difference for 
the subsequent considerations we will restrict to the two-scale case. 

From (|6.4|) we see that for any e > 0, and given hi and /12 there exists D > such that the 
linear combination 

(^D)-"/2( Y e-l^-^il'/"?^+ Y ac;2e-l"-^^l'/"'^) (6.5) 
giGGi g2eG2 

with = Ok for 52 = /iini-|-/i2k, m G Z2, k G S, approximates the constant function 1 with an 
error less than e/2. The idea of constructing the 0-function ()6.1() is to choose for each gi G Gi 
and 52 £ G2 finite sets of nodes ^(ffi) C Ji and £(32) C J2, respectively, and to determine 
polynomials V^^g^ such that 



Y Pj,,,(^-^)e-l^-^jl'/^'^ approximate e'l'^-^'^l'/^'^ , ^ = 1, 2 . 



If the Loo-error of the sums over gi can be controlled, then we get 



with the polynomials 

E (6.6) 

9ieG(xj) 

and 

E a,,e-l'^-s^l'/'^i^ X E vJ^^^)e-\^~^-\'/^''' 
with the polynomials 

Vk= E 092^k,s2 , (6.7) 

92eG(xk) 

where we denote G(xj) = {g : £ ^{d)}- Note that we have to choose the subsets ^(g^) such 
that the sets G(xj) C Gg are finite and nonempty for any node xj € Jg. Additionally, one has 
to choose these sets such that for some ki > and any gi € Gg the ball B{gi, Kihi) contains at 
least one node xj G J^. This is always possible, since Conditions 12.11 resp. 16. II are valid. 

The proposed construction method of does not require solving a large algebraic system. 
Instead, to obtain the local representation of one has to solve a small number of approximation 
problems, which are reduced in the next sections to linear systems of moderate size. 

After this preparation we write G as 

e(x)=(7rD)-"/2 E e-l^-^il'/'^?^ + E 

+ {.D)-i^ E + E 

52eG2 g2eG2 



where 



(y) = (vrD)-"/2{ E Pj,,,(^^)e-|y-yjl'/^ -e-ly-^^/'^^l'/^} (6.8) 



with yj = Xj//if, Xj G Jg. Hence for sufficiently large D 

\ <i.„. ii]„. i 

-h2 



|e(x)-i|<| + E h9i(:^)|+ E ~^92\^92{^) ■ (6-9) 

giGGi 92eG2 



6.2 Construction of Polynomials 

Let us introduce 



(y) := (vri?)-"/2{ ^ p. (^^1 ^eH^"^^!^/^ - e"!^!^/^}, (6.10) 



where S is some finite point set in M". We will describe a method for constructing polynomials 
such that e'^l^l l'-^(y)l foi^ some p > becomes small. In what follows we use the representation 

h 
1-31=0 



Hence by Lemma l5. II 



^ |/3|=0 

and uj can be written as 

..(y) = (vri?)-"/^ ( J2 ^i,f^M^dy) e'l^-^^l'/^ - e-l^l^/^) . (6.11) 

yj6S|/3|=0 

To estimate the Loo-norm of uj we represent this function as convolution. 
Lemma 6.1 Let V{t) be a polynomial and let < Dq < D. Then 

V{d^) e-l«-l'/^ = ci e-l«l'/(^-^") * V{d^) e-l«-l'/^o , 
where * stands for the convolution operator and 

D \n/2 



.7tDo{D-Do) 

Proof. From 



we obtain 



P(5^)e-l«-^l'/^ = p(-ax)e-l«-^l'/^ = ci y"e-l«-*l'/(^-^«)p(at)e-l*-^l'/^« dt . ■ 

Using Lemma l6. II and 1)6. 11() we write uj as 
c.(y) = -^/e-l-*l^/(^-^o)(^ Ec,^5M^/^9t)e-l--.lV^o_e-|t|Vz.o),t (6.12) 

and, by Cauchy's inequahty, we obtain 

l^oo ^ -2 II E E cj,/35MV^50e-l--^lV^o _ e-|t|V^o , (6.13) 



yjGS|/3|=0 

where 

C2 = {7rDor^/\27T{D-Do))-^/'. 
If we define polynomials by 

r^(x) = e l^l'/^05^(\/;D9x)e-l''l'/^« , (6.14) 

then 



yjes 1/91=0 

An estimate for the sum of IwoJ can be derived from 



L2 



D-Do 

Lemma 6.2 Let < Dq < D and denote p = -— — -^^ . Then the estimate 

[D - DqY + DDq 



sup|L^(y)|e''lyl' ^caVoR (6.15) 



is valid, where for c = the quadratic form Q{c) is defined by 

Li 



g(c) = L2iD-Do)\t\yDDo(^^-\t\yDo J2 Cj,^T^(t -yj)e-|t-yjl'/^«)'dt (6.16) 

---^1^1=0 



yj62|/3|=0 

and 

C3 



{27T^Do{D - Do){{D - DoY + DDq)Y/^ 
Proof. Starting with (jlTT^ . using (|6.14|) and 



|x-t|2 



\/^x - V (1 - a)|xp + - — 3:|t|2 
\/a a 



for a > 0, we derive the representation 

^fvl ^1 c-(l-a)lylV(g-go) / -|t-ax|2/a(i?-Do) (l-a)|t|Va(D-Do) 

(7rZ))"/2 y 

X ( E E -SM^ - yj)e-l*-^^l^/^" - e-l*l^/^°) dt . 
yjGS|/3|=0 

Then Cauchy's inequahty leads to 
^(y)e(i-«)MV(^-i5o) 

<C3( L'(^~-W'MD^Do)(^ ^ ^ Cj,^r^(t-yj)e-l*-yjl'/^o-e-l*l'/^«)'dt)'^' ^^'^^^ 
j/„ yjes 1,31=0 

with 

n/4 



C3 = {ttDq 



-n/2 



.27r(L'-L'o)' 
If we choose the parameter a such that 

(l-a)|t|2 |t|2 |t|2 . DDq 



I.e. a 



aiD-Do) Do D ' • ■ {D-Doy + DDo' 

then the right hand side of 1)6. 17() takes the form ()6.16|1 . ■ 

Next we estimate 

r:=minQ(c). (6.18) 

c 

Using ()6.14p . after elementary calculations one obtains 

Q(c) = / ^^iD-Do)\t\VDDo^ ^ J2 Cj,^5^(V^9t)e-l*-''^l'/^« -e-l*l'/^«))'dt 

yjGS |/3|=0 

= (^) ^ (1-2 E E Cj,/3C;3,o(yj,o) + E E E Cj,/3Ck,^c^,^(yj,yk) 

yjeS 1^1=0 yj,ykeS|/3|=0|')'|=0 



with 



The minimum of Q{c) is attained by the solution c = {cj f^} of the hnear system 
^ ^ (^A-clyj ' yk)cj,/3 = Co,-y{0, yk) , Yk G S , < I7I < Lk . 

yj6S|/3|=0 

Then by Lemma 16.21 the sum 



i^ y-_yj V-iy~y,p/D 

D 



approximates e"'^'^/^ with 



yjeE 1^1=0 



y ~_yj \„-|y-yilVg 



<C3e-^|ylVi/2 



(6.19) 



(6.20) 



In the next section we show that (|6.19|) has a unique solution and give an estimate of r. 



6.3 Existence and estimates 

Let us give another representation of the quadratic form Q{c) defined by H6.16|) . Introduce the 
transformed points 

D 



tj = -FT^i ' yj e s 



then, because of 



D-D, 



D-D, 



t t-y; = t-tir + 



I |2 

'yjl 



(5(c) can be written as 



Qic) 



,(D-Do)|yjlV^5g 



E Cj,^r^(t-yj)eH*-*jl'/^) dt. (6.21 



|/3|=0 



Since are polynomials of degree /3, the minimum problem for Q{c) is equivalent to finding 
the best L2-approximation 



mm 



yjeS|/3|=o 



Lemma 6.3 Lei {xj} a finite collection of nodes. For all > the polynomials "Pj of degree 
Li, which minimize 



L2 



(6.22) 



are uniquely determined. 



Proof. The application of Lemma 15.11 gives for 'Pj(x) = c^j^x^ 

1,31=0 



-1-1 = 



7r\"/2 
2 



j l/3|=0 j,k |^|,|tI=0 



where we use the notation 



fi;3,-y(x, y) = Sf3{-d^) S-,{-dy) e-l'^-yl'/^ . 

The coefficients {cj^^} are chosen to minimize ()6.22() . that is the vector {cj^^} is a solution of 
the linear system 



Cj,/3%7(xj,Xk) = eo,-7(0,Xk). 

j \f3\=0 

To show that the matrix of this system is positive definite we use the representation 



(6.23) 



which implies 



%-y(x,y) = (27r)-"/2y" Sf3{-it)S-,i-it) 



g~|t|V2gi(t,x)g-^(t,y) _ 



Let {vj^is} be a constant vector and consider the sesquilinear form 



j-k |/3|,|7l=0 



J.k |/3|,|t|=0 



(2v.)-"/2|eH*lV2| ^ ,.^5M-zt)e^(*-^) 



dt>0. 



j l/3|=0 



The change of integration and summation is valid because the integrand is absolutely integrable 
and the sums are finite. We have to show that the inequality is strict when {fj,/3} 7^ 0. This is 
equivalent to show that 



j l/3|=0 



identically only if all components Vj^f3 = for all j and (3. To this end similar to jl9l Lemma 3.1] 
we introduce the function 



R" j 1^1=0 ]gn 



j l/3|=0 

Li 



j 1^1=0 



Let us fix k and consider /£(x) for |x — x^l < e for sufficiently small e > 0. We have 
5^(5x)e-l''-''jl'/^' ^ as e ^ and 



|/3|=0 1^1=0 

Because of /e(x) = for all e > there exist Eq such that 



t=(x-Xk)/e 



E ^^k,/3 5;3((?x)e-l''-'"'l'/^' = (6.24) 

1/31=0 

for e < eo- On the other hand, 

Li, 

E ^k,^5^(5x)e-l^-^''l'/^' =e-2^''e-l---''lV-^n2L,(e), 

|/3|=0 

where n2L^(e) is a polynomial of degree 2Lk in e with coefficients depending on |x — Xk|. 
Therefore 1)6. 24|) holds for any e > 0, in particular 

Li, 

E ^k,/3 5^(9x)e-l''-^''l' = 0. 

|/3|=0 

Since by 

5/3(ax)e-l^-^>^l' = (x - xk)^e"l^-^>^l' , 
we conclude Uk,^ = for all /9. ■ 

Let now for given S and degrees Lj the coefficient vector c = {cj^^} be a unique solution of 
the linear system 1)6. 19p . To estimate r = Q{c) we denote by G S the point closest to and 
by the degree of the polynomial "P^. 

Lemma 6.4 The minimal value of (|6.18() can be estimated by 

'■jT\n/2 JjL^+l+n/2 |y^|2(L^+l) 



r < 



Proof. It follows from the representation (|6.21|) that 

'n yjeE |/3|=o 
< min [ fp(t)e-l*-*-lV^-e-|tlV^f dt 



2 / PeHi 



with = Dj^/Dq, Zfj, = V Dy ^/ Dq, and 11^^ denotes the set of polynomials of degree L^. 
The minimum is attained when 



1 

V{t) = , y aBHd^) 



|/3|=0 

with the coefficients 



a/3 



J ^ a/2I/3|/3!W2 J ^ ' 



/3!W2 J ' V2l/3|/3!7r"/2 

Integrating by parts, we obtain 



which together with 



leads to 



6.4 Approximate partition of unity with Gaussians 

Now we are in position to prove the main result of this section. Suppose that the nodes {xj} 
are as described in subsection 16.11 and let Gi U G2 be the associated piecewise uniform grid with 
stepsizes hi and /i2- Assign to each grid point S Gi, I =1,2, a, finite set of nodes T,{g(), fix a 
common degree L for all polynomials "Pj in (|6.1|1 and solve the linear system 

E EC^.(^.2^)aM«)^C0.(0,-^) (.25) 

XjeS(9^)|/3|=0 
for all Xk G "^{gi) and < I'/l < L with 

C/3,^(x,y) =5/3(-V^5x)5^(-y^5y)e(^-^o)(l-l'+|y|')/^oe-^l-y|'/2^o , 



Dq < D is some arbitrary positive number. Following (|6.2U() define the polynomials 

L 

/X — X: \ P 

Xj e Ji, 



/X — Xj 



giGG{xj) |/3|=0 '^1^^ 



(6.26) 



926G(xk) |/3|=0 



Theorem 6.1 Under Conditions 12.11 and 16. II on i/ie scattered nodes {xj} /or any e > t/iere 
exisi D > and L such that the function (|6.1|) is an approximate partition of unity satisfying 

|e(x)-l|<e /oraZ/xGR", 

if the polynomials {Vj} of degree L are generated via (|6.26j) by the solutions {c^^/sige)} of the 
linear systems (|6.25|) for all £ Gi U G2- 



Proof. From ()6.9|) we have to show that 

sup( + E 

916G1 g2eG2 
if L is sufficiently large. We start with estimating the first sum 



^92 



LJ 



''^h2' 



e 

< - 
- 2 



(6.27) 



3ieGi 



where 51 = /iim, m G Zi c Z". Using (|6.1U() we can write 



where the points yj in (|6.1Uj) are given by yj = Xj//ii — m, Xj G T,{gi). By Lemmas 16.21 and IH7 

we have 



E h».(^) 



9ieGi 



vr\"/4 ^ _,2 D(i'M^+l+"/2)/2 



meZi 



■p|x/hi-mp 



m 



where x^^ G ^(511) is the node closest to gi = him and is the degree of the polynomial 
Since |x^^ — /iim| < Kihi by Condition 12 . 1 1 and Lf^^ = L for all we conclude that 



E h.^(^ 

9ieGi 



<C3 



2^ Z)o^+^y(LT 



'1 \ " 

sup > 

1)! R" 4^ 



^-p|x//ii-m|2 



(6.28) 



From 



^ (D - Do)2 + L'A 
we see, that for fixed D and Dq 



G (0, D) for any fixed L>o e (0, L>) , 



E 

9ieGi 



if L ^ 00 



(6.29) 



We turn to 



E ^92(1^) 



92SG2 



with g2 = him + /12k, m E k S S. Using H6.10() we have 



,x — uih 



and the points yj in ()6.1Up are given by yj = (xj — m/ii)//i2 — k with xj G Ti{g2). Hence 



E 

92GG2 



<C3 



"92 



TT 



'ii 



2^akk^( 7- k) 



mGZa keS 



n/4 



EE 

meZ2 kG5 



Oke 



-p|(x-m/ii)//i2-k|2 



_f)(LMk+l+"/2)/2 



x^^ - m/ii _ 



/l2 



Here x^^^ G ^((72) is the node closest to §2 = /iim + /12k and L^^^ is the degree of the polynomial 
Vfj,^^g2- By Condition 16. II for fixed D and Dq 



j^{L+l+n/2)/2 



Xp^ - mhi 



ho 



L+l 



< 6{L) ^ if L^oo 



uniformly for all 52 G G2. Hence we obtain 



E' 

326G2 



^92 



UJ. 



''^h2 



<Ci5{L) Y^a^^-p\{^~^h^)/h2--k\^ 
meZ2 keS 



(6.30) 



because of L^^^ = L for all /x^- The sum 

-p|(x-m/ii)//i2-k|2 



Oke 

keZ" 



keZ" 



can be easily estimated by using equation ()6.2() . Setting 

(/i? - /ii)Z) = hjDi - hl/p 

we derive 



Di=D + ^(--D] =D + H' 



h\\p 



D-Dn 



and after some algebra 



^-h^,\k\y{hi-hi)D g-p|x-h2kiv/«i 



rDfl „/2 



VW'/^?^^fl + 0(e 



-w^D^(l-H^)/Di 



Therefore we obtain 

sup V Vake-''l(''-'"'^i)/"^-'^l' <C2Sup V e"!— < C3 

meZ2 ke5 mgZ2 

with some constant C3 depending on D, Dq and the space dimension n. Now ()6.27p follows 
immediately from ()6.29|) and H6.3U() . ■ 

Remark 6.1 It can be seen from 1)6. 28p that in principle the parameter Dq can be any value of 
the interval (0, D) not too close the its end points. In numerical experiments we have not seen 
any significant dependence on this parameter. The choice Dq = D/2 might be advantageous 
because the differential expressions Cf3^-y{x, y) simplify to 



6.5 Numerical Experiments 



We have tested the construction (|6.26l I6.25|) for a quasi-uniform distribution of nodes on M with 
the parameters D = 2, Dq = 3/2, h = 1, ki = 1/2. To see the dependence of the approximation 
error from the number of nodes in 5](m), m S Z, and the degree of polynomials we provide 
graphs of the difference to 1 for the following cases : 

- T.{m) consists of 1 point, L = 1,2 (Fig. ^ and L = 3,4 (Fig. EJ; 

- S(m) consists of 3 points, L = 1, 2 (Fig. © and L = 3, 4 (Fig. I7|); 

- S(m) consists of 5 points, L = 1, 2 (Fig. |HI) and L = 3, 4 (Fig. 0. 

As expected, the approximation becomes better with increasing degree L and more points in the 
subsets S(m). The use of only one node in S(m) reduces the approximation error by a factor 
10^^ if L increases by 1. The cases of 3 and 5 points indicate, that enlarging the degree L of 
the polynomials by 1 gives a factor 10~^ for the approximation error. 

One should notice, that the plotted total error consists of two parts. Using (|6.261 16.25|) we 
approximate the 0-function 

oo 

(27r)-i/2 ^ e-(^-'-)'/2 = 1 + 2 e-'^"^" cos IttJx . (6.31) 

meZ j=l 

Hence, the plotted total error is the sum of the difference between 1)6. 1|) and ()6.31|) and the 
function 

oo 

2^e-^'''^'cos27rjx, (6.32) 
i=i 

which is the saturation term obtained on the uniform grid. The error plots in Figure [7| for L = 4 
and in Figure 1^1 show that the total error is already majorized by H6.32() . which is shown by 
dashed lines. 
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Figure 4: The graph of 0(x) — 1 when S(m) consists of 1 point, L = 1 (on the left) and L = 2 (on the 
right). 



-4*10 



5*10"= 




-1 


-0 . 5 




\) . 5 


1 






-5*10"= 







Figure 5: The graph of 0(x) — 1 when I](m) consists of 1 point, L — 3 (on the left) and L — 4 (on the 
right). 




Figure 6: The graph of 6(x) — 1 when I](m) consists of 3 points, L = 1 (on the left) and L = 2 (on the 
right). 
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Figure 7: The graph of 8(x) — 1 when I](m) consists of 3 points, L = 3 (on the left) and L — A (on the 
right). The saturation term (|6.32() is depicted by dashed lines. 



Figure 8: The graph of 8(x) — 1 when S(m) consists of 5 points, L = 1 (on the left) and L = 2 (on the 
right). 




Figure 9: The graph of 8(x) — 1 when I](m) consists of 5 points, L — 3 (on the left) and L = 4 (on the 
right). The saturation term (|6.32|) is depicted by dashed lines. 
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